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Abstract 

Biophysicists and structural biologists increasingly acknowledge the role played by the 
mechanical properties of macromolecules as a critical element in many biological processes. 
This change has been brought about, in part, by the advent of single molecule biophysics 
techniques that have made it possible to exert piconewton forces on key macromolecules 
and observe their deformations at nanometer length scales, as well as to observe the me- 
chanical action of macromolecules such as molecular motors. This has opened up immense 
possibilities for a new generation of mechanical investigations that will respond to such 
measurements in an attempt to develop a coherent theory for the mechanical behavior of 
macromolecules under conditions where thermal and chemical effects are on an equal footing 
with deterministic forces. This paper presents an application of the principles of mechanics 
to the problem of DNA packaging, one of the key events in the life cycle of bacterial viruses 
with special reference to the nature of the internal forces that are built up during the DNA 
packaging process. 

1 Introduction 

Mechanics has a long and rich tradition of reaching out to other fields. For example, the 
mechanical analysis of thin films has revealed insights into applications ranging from microelec- 
tronics to MEMS (see Freund and Sm'esh (2003)). Mechanics arguments have demonstrated 
that dislocations and surface instabilities play a key role in determining the mechanical and 
electrical properties of such films. While the study of thin films is a relatively nev^^ area to have 
been influenced by mechanics there are others, like geology and metallurgy, that have drawn on 
mechanics for decades. The most recent developments in intersonic crack propagation (Rosakis 
(2002), Rice (2001)) have direct analogues in the study of motion along geological fault lines. 
Geophysicists have relied heavily on data emerging from shock compression of solids to model 
the thermomechanical behavior of planets (see Meyers (1994)). In addition, the impact of dis- 
location mechanics, phase transformations and fracture mechanics on the science of metallurgy 
can scarcely be overestimated. 

Another field that promises to be visited by mechanicians with increasing regularity is 
biology, though the use of mechanical principles for both fluids and solids in understanding 
biological processes is not at all new. For example, slender body theory from fluid mechanics 
has been fruitfully applied to the study of motion in bacteria and other microorganisms (see 
Bray (2001)). In addition, variational principles for determining optimal shapes have been 
used for the study of red-blood cells (see Boal (2002)). Evans and Skalak (1980) show how 
lipid bilayer vesicles can be modeled as shells and membranes with no shear resistance. The 
flagella of bacteria have been idealized as rods in an effort to study how their rotation can 
give rise to propulsion (Goldstein et al. (1998)). The cytoskeletal framework of the cell is 
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studied as a three dimensional network of elements to obtain estimates of the response of cells 
to external forces (see Boal (2002)). In this paper we restrict our attention to problems in 
biological nanomechanics with the observation that there are a range of fascinating connections 
between mechanics and biology at larger scales as well. One of the most intriguiging features 
of problems in biological nanomechanics is that they involve a rich interplay of statistical and 
deterministic forces as will be seen below. 

It is only recently, however, that structural biologists have begun to understand the signifi- 
cance of mechanical properties for the macromolecular processes that sustain life. The role of 
bending stiffness of DNA in gene regulation (Widom (2001)) and the stiffness of substrate in cell 
migration (Lo et al. (2000)) are two examples. Such insights have been garnered with the aid of 
sophisticated experiments that have made it possible to apply piconewton forces on nanometer 
size objects and measure the resulting deformations. The data from these experiments point 
to a rich interplay of thermal and chemical forces with electrical and mechanical forces. In 
this paper we show how experimental insights can be combined with mechanical principles to 
construct a simple and quantitative mechanical theory of DNA packaging in bacterial viruses 
that suggests a new round of experiments. 

2 Mechanics and the Viral Life Cycle 

Viruses have been studied extensively (see Alberts et al. (1997)) in the last three decades not 
just with the goal of understanding their pathogenic nature but also for gaining insights into 
structural biology of proteins and nucleic acids and for studying gene regulation. Even though 
viruses are perhaps the simplest entities from a biological perspective, they are nanotechnologi- 
cal marvels that embody a wealth of physics at the nanometer scale. The principles operating at 
these scales are being rapidly unravelled through ingenious experiments, such as that by Smith 
et al. (2001) on DNA packaging in the (p29 virus, which serves as the primary motivation for 
the theoretical analysis presented below. 

Before going into the details of the experiment it is useful to get a glimpse of the life-cycle of 
a virus (see fig.^. We are especially interested in the class of viruses known as bacteriophage 
and which infect bacteria such as the well-known E.coli. In simplest terms, a bacteriophage 
is nothing but a protein coat, the capsid, filled with nucleic acid such as DNA or RNA. The 
capsids of bacteriophage attach to the surface of the bacterial cell which is under attack. The 
genetic material of the virus is ejected into the bacterium leaving the capsid behind. Though the 
process of DNA ejection is not completely understood, it is argued that a mature bacteriophage 
capsid is highly pressurised (see Smith et al. (2001)) and the genome is forcefully released into 
the host cell whose contents are effectively at a much lower pressure. 

Once the viral DNA is inside the bacterium it hijacks the protein production machinery of 
the bacterium to synthesize its own proteins such as those that will ultimately make up the 
capsid. Copies of the viral DNA are also made. Here again, experiments have shed light on the 
mechanics and kinetics of the processes of transcription of the genome and its translation into 
proteins. The process of transcription is mediated by a large protein called RNA polymerase. 
This protein attaches to the DNA and then moves along its length transcribing the DNA into a 
molecule known as mRNA. As it moves along, the protein exerts forces on the DNA molecule. 
These forces have been measured experimentally (Wang et al. (1998)) and it has been found that 
the protein can be stalled by exerting a force of roughly 25pN. In fact, the rate of transcription 
is heavily dependent on the force exerted. The type of data emerging from such experiments 
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unmistakably points to a plethora of mechanics problems still to be resolved at the level of 
macromolecular assemblies. 

After a sufficient quantity of the viral proteins have been synthesized they begin to self- 
assemble into hollow capsids. In the case of the (j)29 virus, the capsids have a portal at one end 
with an attached protein motor. The motor proteins identify the appropriate end of the as- yet 
unpackaged viral DNA and push it into the capsid. As more of the DNA gets pushed into the 
capsid the motor has to perform work against an increasing resistive force due to confinement 
of the packaged DNA. This causes the motor to slow down as more DNA fills up, but ultimately 
all the DNA is packed and the remaining proteins attach themselves to the capsid thus making 
it ready to infect another bacterium. Once the capsids are fully assembled an enzyme is released 
that breaks up the bacterial cell membrane and releases the mature viruses so that they can 
repeat their evil action elsewhere (see Ptashne (1992)). 



Infection 




Figure 1: The life cycle of a virus. The critical step which is the focus of this paper is the 
packaging of DNA into the assembled capsid. The packaging process involves large forces on 
the order of several tens of pico-newtons and it proceeds at high rates, packing the full length 
of the genome in only a few minutes. 

As noted above, a critical step in the life cycle (see fig. ^) of the virus is the packaging of 
DNA into the capsid by the portal motor. This step was the focus of the experiment by Smith 
et al. (2001) in which they used optical tweezers to measure the force exerted by the portal 
motor of the 4>29 bacteriophage on the DNA as a function of the amount of DNA already packed 
in the capsid. One end of the DNA was attached to a silica bead which was held in an optical 
trap while the motor tugged on the other end. The capsid itself was attached to a second 
silica bead that was held by a glass pipette thus immobilizing the capsid during the packaging 
process. Through a series of careful experiments, the average force exerted by the portal motor 
was measured and plotted against the fraction of packed DNA. The rate at which the DNA is 
packed was also measured as a function of the fraction packed. It was found that during the 
early stages of packing the motor packs at a rate of roughly 100 base pairs per second which 
correponds to 34nm per second. However, the rate drops considerably after about 70% of the 
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DNA has been packed and it slows to a crawl as the packaging nears completion. We note that 
the experiment was conducted in vitro in a solution which contained 50mM NaCl, 5mM MgCl2 
and tris-HCl buffer at pH 7.8. These are not the conditions prevailing in a bacterium and one 
of the interesting outcomes of the type of model being described here is the recognition that the 
ionic strength is a key control parameter in governing the physics of DNA packaging. There are 
a variety of interesting mechanics issues to be confronted relating to those stages of the viral 
life-cycle involving both DNA packaging and ejection. One such issue involves determining 
what is the maximum sustainable pressure in a viral capsid (Purohit et al. (2003)). Another 
concerns determining the amount of DNA ejected by a bacteriophage as a function of the 
osmotic pressure in the exterior solution (Evilevitch et al. (2003)). 

This paper puts forth a simple mechanical model of the packaging process in an endeavour 
to understand the nature of the forces being exerted on the DNA. Previously, we described 
results of an analysis based on a continuum approximation to the mechanical model of the 
packaging force (Purohit et al. (2003)). Here we analyze a discrete version of that same model 
which reveals interesting new effects which could be accessible experimentally. 

3 The Free Energy Function for Packed DNA 

In order to estimate the energetics of viral DNA packing, wc take a structural hint from the 
cryo-electron microscopy images of capsids (see, for example, Cerritelli et al. (1997)) in various 
stages of the packing process. It has been observed for several different viruses that the packed 
DNA is arranged in a series of circular hoops (or a spool). Each hoop finds itself at the center 
of a hexagon formed by nearly parallel segments of the same overall strand. The exception are 
the hoops hugging the surface of the capsid, or those at the inner surfaces of the spool, which 
are surrounded on average by three nearest neighbors. We note that the geometry is more 
subtle than the hoop packing adopted here since the packing involves a helical pitch and hence 
the adjacent strands are probably not perfectly parallel. Following earlier work (see Riemer 
and Bloomfield (1978), Odijk (1998) and Kindt et al. (2001)), the basic idea is to write a free 
energy function which characterizes the energetics of structures like those described above as 
a sum of elastic terms and interaction terms. It is well known that DNA has a considerable 
bending stiffness at length scales of a few tens of nanometers. This is evident from the fact that 
the persistence length of DNA is about 50nm. The persistence length is defined as the distance 
over which the tangents to a fiber are correlated. The more flexible a fiber is the smaller its 
persistence length. Polymer physicists use the persistence length as a measure of the stiffness 
of a chain (see Doi and Edwards (1988)). Since the 

persistence length of DNA is of the same order as the size of viral capsids we expect the 
bending energy to play an important part in the packaging process. Such a description of the 
elasticity of DNA is now quite firmly established in literature and is found to be useful in 
understanding key biological processes related to genetic regulation and DNA condensation in 
chromosomes of eucaryotic cells. 

Beyond its explicit mechanical properties another relevant characteristic of DNA is its strong 
acidity. In particular, in an aqueous solution the backbone of the molecule is highly charged, 
with two units of negative charge per base pair. As a result there arc electrostatic penalties 
involved in bringing two strands in close proximity. In the salty solutions found in cells the 
picture is more complicated because charged ions screen these interactions. Another compli- 
cation arises from the fact that close packing involves the ejection of water molecules from 
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the vicinity of the strands, which make sizable contributions to the free energy. All of these 
complex interactions of the DNA with the ambient solution and with itself result in an effective 
interaction energy. The importance of these effects is immediately evident when we look at 
the physics of DNA condensation in a solutions containing trivalent and tetravalent cations. In 
particular, DNA is known to form hcxagonally packed toroidal structures when immersed in a 
solution containing these ions (Raspaud et al. (1998)). The aim of the calculations to follow is 
to model viral packing, taking into account both elastic and interaction energies. 

3.1 Elastic Contribution to the Free Energy 

We begin by idealizing DNA as an elastic rod capable of sustaining bending deformations. The 
centerline of the rod is parametrized by the reference coordinate x. The bending energy stored 
in length / of this rod is given by 



where R{x) is the radius of curvature at reference position x along the rod and k is a bending 
modulus which we assume to be independent of the reference position x though we note that 
the 

sequence dependence of the elasticity of DNA is one of its most 

intriguing features. For a rod bent into a hoop, the radius of curvature R(x) is just the 
radius of the hoop R, and the length is / = 2ttR, and so the integral reduces to 



The modulus k is usually taken to be EI where E is the Young's modulus of the material 
and / is the moment of inertia of the cross-section. However, it is convenient to express the 
elastic properties of molecules differently since thermal oscillations play a dominant role in their 
deformation. In particular, the deformation at a given point on a molecule may be uncorrelated 
with the deformation at another location because of the randomness caused by these thermal 
motions. As noted above, the length at which the importance of bending energy is comparable 
to that due to thermal vibrations is called the 'bend persistence length' and is denoted by 
^p. More precisely, is determined by balancing the bending energy and the thermal energy 
and results in «; = ^pksT, where ^p is the persistence length, T is the temperature and 
is Boltzmann's constant; ^p = 50nm for DNA and ksT = 4.1pNnm at T = 300K. Using this 
description of the elastic properties of DNA and adding up the contributions due to all of the 
hoops in the spool, the total bending energy is given by 



where N{Ri) is the number of hoops in a column of radius Note that we neglect the fact that 
the DNA is packed in a helical arrangement with an associated contribution to the curvature 
due to the helical pitch. As noted in our earlier paper (Purohit et al. (2003)), this effect is 
negligible. For analytical simplicity the discrete expression given above can be converted into 
an integral and written as 




(1) 




(2) 




(3) 
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where ds is the spacing between adjacent hoops, Rout is the radius of the capsid and R is the 
radius of the innermost set of hoops. The factor of appears since the spacing between two 

adjacent columns of hoops is ^dg in a hexagonal array. The integral approximation has the 
virtue that it leads to closed-form analytical expressions for the energies and forces in the viral 
packing problem, at least for simple geometries like the cylinder and the sphere. An interesting 

set of conclusions to be discussed below are the differences between the discrete and continuum 
descriptions of this problem, and their possible experimental consequences. 



3.2 Free Energy of Interaction 

In addition to the role played by elastic bending we must also consider the effect of the inter- 
action between adjacent strands of DNA. For this we appeal to the experiments of Parsegian 
et al. (1986) and Rau et al. (1984). In these experiments, osmotic pressure was applied on 
hexagonally packed DNA and the interstrand separation was measured as a function of the 
pressure for a variety of different solvent conditions. For conditions comparable to those of 
Smith et al. (2001), it was found that the osmotic pressure could be related to the interstrand 
spacing as 

p(4) = Foexp(--^), (5) 

where Fq is a constant whose magnitude depends on the type and strength of the ionic solution, 
and c = 0.27nm is a decay length which is roughly constant over a range of ionic conditions. 
For a solution containing 500mM NaCl at 298K, measurements reveal Fq = 55000pN/nm^. To 
this piece of empirical evidence we add the assumption that parallel strands interact through a 
pair potential v{ds) per unit length and that interactions are limited only to nearest neighbors. 
Given this description of the total energy and the measured pressure vs interstrand spacing, we 
can deduce the functional form of v{ds). 

In particular, consider N parallel strands of length I each packed in a hexagonal array with 
a spacing dg. The total interaction energy of this arrangement is 

E = SNlv{ds), (6) 

where the factor of 3 appears because each strand interacts with 6 nearest neighbours (ignoring 
surface effects) and we multiply by 1/2 to account for double counting. The total volume of 
the assemblage is obtained by adding together the (prismatic) volume occupied by each strand, 
thus 

V = N^dll. (7) 

We now use the thermodynamic identity p — — to relate the expression for pressure obtained 
from experiments to our simple model. We note that dE = SNl-§^dds and dV = VSdgNldds, 
so that 

«*' = -^ = ;^*>*- 

where f{ds) is the force per unit length on the strands. We now substitute the experimental 
result, namely p{ds) = Fq exp(— ^) and solve the differential equation above with the boundary 
condition v{oo) = 0. This results in the potential 

v{ds) = -^Foic^ + cds) eM--)- (9) 
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If this result is now exploited in the context of eqn. (jH)) the interaction energy can be written 
as 

Eint = VSFoic^ + c4)iexp(--f ), (10) 

where L = Nl is the total length of the strands. 

Now we are in a position to write the total energy of the packaged DNA as a sum of the 
bending energy and the interaction energy as follows 

E{L, ds) = V3FoL{c^ + cds) eM--) + ^pkBT^,^^^- (H) 

I 

In the continuum approximation this takes the form 

E{L,ds) = V3FoL{c' + cds)eM--) + ^^^^^ T"" ^dr. (12) 

We still need to express the length L in terms of the radii of the hoops. To do so, we simply 
add up the accumulated length of all the hoops as, 

L = ^2-KR^N{Ri). (13) 

i 

In the continuum approximation this can be rewritten as 

L = / rN{r)dr. (14) 



Vsds Jr 

With both the elastic and interaction contribution to the total energy of packed DNA in hand, 
we turn now to a concrete investigation of the implications of this model. 



4 Geometry and Energetics of Packed DNA 

The calculation of the free energy associated with the DNA packed in a partially filled capsid 
is predicated on key insights gained from experiments. We assume throughout that the DNA 
adopts an inverse-spool geometry like that described above, with the proviso that the class of 
minimum free energy structures we consider is constrained to the spool-type and that there 
are perhaps more complex structures with lower free energy. We also assume that the strands 
have local six-fold coordination. As a result, the only variable that remains to be determined is 
the spacing ds between the strands as a function of the amount of DNA packed. In the regime 
where the interaction between the DNA strands is repulsive the interaction term dictates that 
the strands be as far apart as possible. On the other hand, since the outer dimensions of the 
capsid are fixed, larger spacings imply smaller radii of curvature which results in a steep rise in 
the bending energy cost. The competition between these two contributions to the free energy 
determines the geometry of the packed DNA (Reimer and Bloomfield (1978), Odijk (1998), 
Kindt et al. (2001)). We illustrate this through two examples and also make a comparison with 
the experimental data of Smith et al. (2001). 

For simplicity we begin by considering a cylindrical capsid which is a first approximation 
to the geometry of the 4)29 virus considered in the experiments of Smith et al. (2001). In 
reality, the 029 virus is shaped like a hollow oblong spheroid (Tao et al. (1998) and Wikoff and 
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Johnson (1999)) which has an outer radius of 21nm and a height of about 54nm. Its protein 
coat is about 1.6nm thick (on average) so that the radius of the capsid volume available to the 
DNA is Rout = 19.4nm while its height is roughly Sl.Onm. The capsid has other interesting 
geometrical features, but for the purposes of the present calculation we idealize it as a cylinder. 
The dimensions of this idealized cylinder are given in fig. |21 where the key point is selecting the 
height of the cylinder so as to gaurantee that the volume available for packing is identical to 
that of the real virus. The figure shows a second geometry that will be used later as a more 
refined model of the (j)29 capsid. 



Ro„, = 1 9.40nm R„„, = 1 9.40nm 

h = 37.87nm z=12.0nm 




(a) Capsid in 3D (b) Cylinder (c) Capped cylinder 



Figure 2: Model representations of viral capsids. (a) Capsid modeled as a cylinder with hemi- 
spherical caps shown in perspective. The curved rods inside represent the packed DNA. (b) The 
cylindrical geometry shown in cross-section with height adjusted to correspond to the correct 
internal volume of the capsid. The gray circles represent the hexagonally packed DNA hoops, 
(c) Idealization of (p29 as a cylinder with hemispherical caps shown in cross-section. 

The (p29 capsid encloses about 6.6;um of double-stranded DNA corresponding to a genome of 
19.3 kilobase pairs. As noted above for our cylindrical model we take the height of the cylinder 
as z = 37.87nm to ensure that the internal volume is the same as that of the actual capsid. 
For such a capsid, the number of hoops at radius Ri is given by N{Ri) = h/ds where dg is the 
separation between the strands and h is the height of the cylinder as shown in fig. |21 One of the 
key themes of the present paper is the comparison of the discrete and continuous treatments of 
the structure and energetics of viral DNA packing. To that end, it is of interest to examine the 
length of packaged DNA in the discrete setting. Assuming that there are n columns of stacked 
hoops (see fig. |21 where we show two completed columns of hoops, i.e. n = 2), the length of the 
packed DNA is 

L = 2vr— (^Rout + Rout — ~^'^s + Rout — + ••■ + ^out — (n — l)-^^^^ 

This equation results from adding up the accumulated length of the hoops at each radius as 
demanded by eqn. (|13|) . The hexagonal coordination of the strands appears once again through 
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the distance between adjacent columns of hoops which is ^dg- As a concrete example, we use 
the numbers for the (j)29 phage and compute how many columns of stacked hoops it contains 
when fully packed. X-ray measurements on packed i?i29 indicate that the spacing between the 
strands in a fully packed (j)29 capsid is about 2.8nm (Earnshaw and Casjens (1980)). The length 
of the genome is L = 6584 nm, and eqn. (|15|) implies n = 5.61 corresponding to 5 columns of 
completely stacked hoops and one innermost column of hoops that is only about 60% full. There 
is a hollow cylindrical region inside whose radius is Rm = 7.28nm. This calculation reveals that 
there are not a large number of columns of hoops even in a completely filled capsid and as a 
result, caution should be exercised in making the continuum approximation in converting the 
discrete sums in equations @ and ()14() to integrals. 

We now turn to the more detailed aspects of the geometry and see how by minimizing the 
free energy we can obtain experimentally falsifiable insights into the packing process. 

4.1 Interstrand Spacing in Packed DNA Capsids: Discrete Model 

We mentioned earlier that the spacing between the strands is determined by the competition 
between bending and interaction energies. We demonstrate this minimization procedure with 
a very simple calculation. The model as set forth above has one free parameter, namely, Fq^ 
the parameter characterizing the strength of the repulsion between adjacent DNA strands. We 
fix Fq = 225500pN/nm^ since, as shown below, this leads to the best fit to the experiment of 
Smith et al. (2001). We note that the ion concentration in their experiment is small compared 
to the values reported by Rau et al. (1984), which implies that the repulsive forces will be 
stronger than those measured by Rau et al. (1984) resulting, in turn, in a larger value of Fq. 
We continue with the cylindrical capsid geometry with the same internal volume as the i;^29 
phage. The goal is to determine the spacing dg between the strands for a given length L of 
genome packed. To that end, we minimize the free energy with respect to ds at fixed length L. 
In particular, we set ds to a fixed value in eqn. (|15() and then solve for the number of columns 
of hoops this implies while holding the length L fixed. We then have all the information needed 
to calculate the energy using eqn. ()11|). We repeat these steps for a range of different values of 
ds to get a plot of energy vs dg. Finally, we use this plot to determine the dg corresponding to 
minimum energy at a given length of DNA. One such plot is shown in fig. |31 The interaction 
energy decreases as the spacing increases. This is expected since the interaction energy is a 
decaying exponential and increasing dg implies that the interactions come from deeper in the 
tail of the pair potential. The bending energy rises as the spacing increases since the radii 
of the innermost strands become successively smaller. The total energy is dominated by the 
interaction part at smaller spacings and by the bending energy at higher spacings. For the 
case shown in fig. |21we can clearly identify a minimum corresponding to dg = 3.81nm in the 
total energy. More generally to find the trends for all values of the length packed, this same 
procedure is repeated and then we plot the optimal dg (at minimum energy) against the fraction 
of DNA packed for the (j)29 phage as shown in fig. |l]for three different values of Fq. We note 
that in the later stages of packing all the curves collapse onto a single curve, though we also 
note that the discrete and the continuous curves have different values in the large packing limit. 
As shown in our earlier work (Purohit et al. (2003)), dg in the large packing limit is dictated 
by packing the DNA in such a way as to maximize the distance between adjacent strands and 
should be seen as a geometric limit. 

On the basis of the simple model presented above, fig.|^provides several distinct predictions. 
First, on a gross level, it gives the history of the variation of interstrand spacing during the 
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Figure 3: Energy vs interstrand spacing for L = 4.0/xm in a cylindrical capsid with Rout = 
19.4nm and z = 37.87nm and with repulsive parameters Fq = 225500pN/nm^ and c = 0.27nm. 
A minimum is evident at dg = 3.81nm. The sharp turn in the bending energy at dg = 3.81nm 
is the point where a new column of hoops with smaller radius begins to form. 

packing process. This is a verifiable prediction since it is possible to perform the viral packaging 
reaction for different ionic concentrations of the ambient solution and for different values of the 
genome length. Such experiments have been carried out for the T7 phage for three different 
genome lengths by Cerritelli et al. (1997). Similar experiments have also been performed for the 
A-phage (see Earnshaw and Harrison (1977)). The second more subtle outcome of the model is 
the possibility that there are actually discrete effects in the packing process due to the packing 
of the DNA at a finite set of discrete radii. It would be of interest to determine whether these 
effects in the packing spacing (and a related force signature to be described below) are present 
in experiments. 

4.2 Interstrand Spacing in Packed Capsids: Continuum Model 

One of the interesting features of the calculations presented here is the ability to contrast 
discrete and continuous descriptions of the energetics of DNA packing. To that end, we now 
repeat the energy minimization argument that was used to determine the interstrand separation 
dg in the discrete setting in the continuum approximation. In particular, we minimize the free 
energy given in eqn. (|T2|) with respect to dg, but subject to the constraint that L is constant. 
In concrete terms this amounts to 

- = -4FoLexp(--)-^=^y^ ^'^+7!^ A ^''^ 
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Figure 4: Spacing between the DNA strands as a function of the amount of DNA packed for 
a model cyhndrical capsid. The spacing between strands in the early stages is larger when 
the magnitude of the interaction energy is larger. In the later stages of packing the spacing 
is dictated by geometrical constraints. The points correspond to the discrete calculation for 
different ionic conditions and the full curve corresponds to the continuum calculation. 



where we have used the constraint 4t- = which takes the form 

dda 

dR _ 2 j^°-'R'z{R')dR' 

dds ds Rz{R) ■ ^ ' 

We substitute for L using eqn. (|14() and then set ^ = to get the following general equation 

/^P , ds. CpknT CpkBTji;-'ziR')/R'dR' 

v3Foexp( ) = -77775 ^ % • (18) 

" R^dl dl j^"-' R'z{R')dR' 

This equation represents a competition between the interaction terms (on the left) and the 
bending terms (on the right). Note that the effect of the capsid geometry appears through the 
second term on the right hand side of eqn. (|18j) . For a cylinder, z{R!) = h, a constant, and 
eqn. H18() takes the form 

/-^ , ds, ipkBT 2CpkBT log(%i) 
V3Foexp(— ) = ^T^T^o ^ 



c R^d^ dl Rlut-R^ 



(19) 



Eqn. H14|) is similarly specialized for a cylinder to yield 



^=^2^^out-R')- (20) 
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Eqns. (fTO|l and (I^U)) feature two unknowns, R and dg, and we solve for them using the Newton- 
Raphson method. The result of this computation is a history of the interstrand spacing ds as 
a function of the length L of DNA packed. This has been plotted as a thick line in fig. 21 for 
Fq = 137500pN/nm^. It is evident that while the continuum version captures the trend quite 
well, it underestimates the value of dg for the cylindrical geometry in comparison with the more 
exact discrete version. Also, details captured by the discrete version, such as the series of steps 
and jumps, are smoothed out in the continuum approximation. 



5 Forces during DNA packaging 

One of the key outcomes of the model presented here is the internal force that builds up 
during DNA packing. The force vs percent packed curves should be seen as a second set 
of experimental implications for this mechanical model which complements predictions of the 
spacing dg vs percent packed. Indeed, once we have the spacing dg as a function of the amount 
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Figure 5: Variation in the stored energy as a function of the length of DNA packed for a model 
cylindrical capsid. The total energy rises as a function of the amount of DNA packed and with 
increasing value of Fq. Also shown as a thick line is the curve obtained from the continuum 
theory. Because the continuum version of dg is smaller than that obtained from the discrete 
model, there is a corresponding increase in the interaction energy which causes the continuum 
energy to be larger than its discrete counterpart. 

of DNA packed we can determine the energy using eqn. (|11|) or (|12|) depending on whether we 
have adopted the discrete or continuum version of the free energy. The variation of the free 
energy has been plotted for the three Fq values in fig. El Recall that the choice of Fq reflects 
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the ionic conditions present during the packaging process. It is evident that the total energy- 
increases at any given value of the length packed as the strength of the repulsive interactions 
(the value of Fq) is increased. Figure also shows the energy obtained using the continuum 
approximation which overestimates the energy but correctly captures the trend. The force 
required to pack DNA is simply the derivative of the energy with respect to the length packed. 
We have calculated this derivative numerically and the results are plotted in fig. |H1 

The behavior of the force curves is very much as expected, with higher forces for larger 
values of Fq. We note that the strong dependence of the maximum packing force on Fq is 
amenable to experimental observation and is worth further attention. In particular, we note 
that the maximum packing force can change by as much as a factor of two depending upon the 
ionic conditions (and the related value of Fq). Perhaps the most striking feature of the force 
curves is the appearance of steps. The same steps are seen in the plots of dg as a function of 
fraction packed. They occur when a new column (stack) of hoops starts to form on the inside 
of the spool. Curiously enough, the experimental force curve reproduced in fig. [3 also seems 
to show discontinuities of the slope at 4-5 different points, which might correspond to the five 
columns of completely stacked hoops needed to fill the (/>29 capsid with its DNA. However, to 
establish that this actually occurs during the packaging of a virus would require a more detailed 
analysis of the kinetics of packing and more careful experiments. The presence of such steps 
in the model results from implicitly assuming that the dynamics of packing follows the sort 
of organized column-by-column hoop packing described above. Therefore, the observation of 
such steps in experiment would be not only a confirmation of the model, but a provocative hint 
about the dynamics of packing. 



6 Application to the 029 virus 

As a final calculation we go beyond the analytic simplicity offered by the cylindrical capsid 
and consider the second approximate geometry shown in fig. |21 namely, the capped cylinder. 
Though the cylindrical geometry offers considerable insight into the mechanics of viral packaging 
it clearly represents an oversimplification of the (/>29 capsid. The (p29 capsid is shaped like an 
oblong spheroid and we approximate it as a cylinder with hemispherical caps in order to make 
comparisons with the experimental results of Smith et al. (2001). The radius of the cylinder 
(and the hemispheres) is taken to be Rout = 19.4nm (see fig. The height of the cylindrical 
portion (also called the 'waist') \s z = 12.0nm. It is useful to study this geometry because one 
can obtain expressions for spherical capsids (which are very good approximations to icosahedral 
viruses) merely by setting the height of the waist to zero. The expression for N{Ri) is now a 
bit more complicated and is given by 



z + 2\ R^ f — Rf 
NiR,) = ^ -. (21) 

Using this formula in eqn. 1)13(1 results in the following expression for the length of DNA packed: 



z + 2JdsiV3Rout - jds) ^3 
L = 27r -Rout + ^ {Rout - \ds) + ... 



z + 2J{n-l)ds{V^Rout-l{n-l)ds) ^ 
+ ^- 2 " " 1)^^^) • (22) 
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Figure 6: Internal force as a function of the amount of DNA packed for a model cylindrical 
capsid. The plot bears the signature of both the total energy (fig. ISJ and the ds (fig. © curve. 
Note that the maximum force (at 100% packing) is roughly proportional to Fq as shown in the 
inset. 



Unlike the cylindrical geometry, this expression is not analytically tractable. However, the same 
procedure can be used for determining the interstrand spacing, and the energy and force as a 
function of percent packed. The results of numerical calculations for this geometry, both in the 
discrete and the continuous setting, are plotted in fig. [Tj We find that Fq = 225500pN/nm^ 
results in a good fit to the experimental data. Moreover, the continuum and discrete versions 
of the theory are in better agreement for this geometry than they were for cylindrical capsids. 

It is of interest to examine the effect of ionic concentration on the forces exerted by the 
portal motor on the DNA. The ionic concentration affects the value of Fq and we have already 
seen that larger values of Fq imply larger values of the force. In the inset to fig. El the maximum 
force as a function of the parameter Fq is shown. This information can be potentially useful 
in the following way. It is known (see Smith et al. (2001)) that the portal motor of the (/)29 
virus stalls at a force of 57pN. In other words, if one were to conduct the packaging reaction 
in a solution with ionic conditions such that Fq > 250, OOOpN/nm^ then the genome would not 
be completely packed since the motor would stall prior to complete packing. We are hopeful 
that experiments carried out with different ionic concentrations would permit an investigation 
of such predictions. We have made similar calculations to those presented here for viruses other 
than 029 and find interesting variations in the maximum packing force from one virus to the 
next which should be similarly accessible experimentally. 
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Figure 7: Force and interstrand spacing as functions of the amount of DNA packed in a capsid 
idealized as a cylinder with hemispherical caps. The hexagons correspond to the experimental 
data of Smith et al. (2001). The thick lines are results of the continuum model and the circles 
connected with the thin line are obtained from the discrete model. A good fit to their data is 
obtained for Fq = 225500 pN/nm^. 



7 Conclusion 

Recent advances in the development of tools for single molecule manipulation have permitted 
the direct mechanical investigation of a variety of biological processes. One intriguing recent 
example involves the direct measurement of mechanical forces during DNA packaging in viruses. 
We have developed a simple, analytically tractable model which responds to these experiments. 
The model emphasizes the role of elastic and electrostatic effects in the packaging process and 
is in good quantitative agreement with experimental data. More importantly, the model makes 
a series of predictions that we hope will inspire and guide future experiments. In particular, we 
predict the appearance of steps in the force versus percent of DNA packed curve. There may 
be some evidence of steps in the existing experimental data, but further experimentation is 
needed. The model also predicts a growth of the maximum force required for packaging as the 
ion concentration is decreased; the underlying cause being the stronger electrostatic repulsion 
between the DNA strands inside the capsid. This points to a new series of experiments that 
would look for incomplete packaging of the viral DNA at lower salt concentrations than used 
previously. An alternative approach would be to work at the same ion concentration but with 
longer DNA constructs. 

The biggest shortcoming of the model as presently stated is that it does not address dynam- 
ical issues associated with the packaging process. The fundamental experimental observation is 
that as the packaging process proceeds the packaging rate falls from the initial 100 base pairs 
per second to zero. We note that the appearance of force steps described in this paper implicitly 
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assumes a precise dynamical pathway resulting in a helical packing of the DNA in a series of 
helices of ever decreasing radius and are hopeful that the future work will shed further light on 
the dynamics of viral packing. 
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